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Abstract 

The propose of this study is to study Physics Learning Achievement, projectile motion, using the Mathematica program 
of Faculty of Science and Technology Phetchabun Rajabhat University students, comparing with Faculty of Science and 
Technology Phetchabun Rajabhat University students who study the projectile motion experiment set. The samples are 
Faculty of Science and Technology Phetchabun Rajabhat University Technology students, studying in the first semester 
of academic year 2011, consisting of 50 samples who study by using the Mathematica program and 50 controlled 
samples, studying by using the projectile motion experiment set. The instrument used for collecting data is the 30 items 
achievement test of the projectile motion. Analyze to compare the achievement between the samples and the controlled 
samples by t-test (Lord, F.M. 1967; Airasian, P.W. and Madaus, G.F., 1972) independent sample. The result finds that 
the achievement of Faculty of Science and Technology Phetchabun Rajabhat University students who study projectile 
motion, using Mathematica is higher than the achievement of Faculty of Science and Technology Phetchabun Rajabhat 
University students who study projectile motion, using projectile motion experiment set at 0.01 level 
Keywords: projectile motion, Mathematica program , pre-test and post-test, t-test, Learning-Classroom 


Introduction 

Nowadays, science and technology play the important role in human living and develop the 
nationwide extremely in teaching management. Physics is the subject for every student who wants 
to pass only, there are very few students who intuitively want to know and understand physics, 
especially the Mathematica (Bauman, 1996) program, Faculty of Science and Technology 
Phetchabun Rajabhat University students don’t pay intention to experiment. So the researcher gets 
interested in using Mathematica program (Bauman, 1996; Bauman, 2005) which is Mathematica 
program (Wolfram, 1994; Maeder, 1991) as the aid for projectile motion (Arya, 1990; Ginsberg, 
1995; Fishbane et al. 2005; Bueche and Hecht, 2006; Marion and Thornton, 1995) experiment to 
assist Faculty of Science and Technology Phetchabun Rajabhat University students to experiment 
the projectile understandingly, to be able to learn and revise it all the time, to learn effectively and 
to solve the students unequal knowledge basic problem. 

In order to study Physics learning achievement, the projectile by using the Mathematica 
program at Faculty of Science and Technology Phetchabun Rajabhat University students. 

The achievement of the Technology students, studying projectile, using the Mathematica 
program is higher than the achievement of the Technology students, studying projectile, using 
projectile motion experiment set. 
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Projectile Motion 

Projectile motion is two - dimensional motion, moving on horizon and vertical 
simultaneously. It is speedy motion because of the steady gravity, being near the surface of the 
earth, whereas there is no speed on the horizon motion owing to lacking of force on horizon, so the 
speed of the horizon movement is steady. The moving route of projectile motion will be parabola 
curve. The projectile motion as the following will be considered without braking power and other 
results such as, rotation, the change of material figures. It is considered only the particle motion, 
which is speedy on only vertical direction. When the motion is specified on one plain, two - 
dimension motion is specified as x-y coordinate system. The origin is zero as x-axis, being on y 
plain. A y-axis is vertical, the gravity acceleration is on -y, or the acceleration component is 


a x = 0 

a y = -y 


( 1 ) 


According to the motion equation, we obtain the velocity on x-axis and y-axis is 


Using equation (2) into 


v ,(0 = Uo 

v y (0 = v yo -gt 


( 2 ) 



x 0 + 



y( t ) = yo + ( v y ( t ) dt 

Jt o 


( 3 ) 


And use the origin to be start point or x 0 = y 0 = 0, we get the particle position for projectile is 

V J 

M 1 2 (4) 

t(0= V" 2 gt 

Restricting t of both equators in (4), the relation between coordinate x and y is (Arya, A, 1990; 
Ginsberg, E.S. 1995) 


g 


y = — =-j— x h—-— x 


2v: 


2 , V y0 

x +■ 


a0 


V 


a0 


( 5 ) 


The relation between x and y of equation (5) is the motion a parabolic trajectory, crossing at x 
axis on 


x 0 = 0 

R 2 U0 V y0 


g 

2 


( 6 ) 
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R Interval, called the range of projectile is the interval that particle move furthest on the plain. 
When the derivative of y in equation (5) comparing with x, the value of derivative is zero to obtain 

the x interval to be —. The interval of y is highest; the value is equal (Fishbane et al. 2005) 


2g 


The figures of parabola mentioned in equation (5) to equation (7) shown as the figure (1) R 
is replaced as the interval on x- axis in the first equation (4), it is the time that the particle being in 
the air. 



Fig. 1. A projectile moving under the force of gravity is at its maximum height when v y = 0 . 

The motion on the plain is equal 


(r - 


2v 


ro 


g 


(0 


— is replaced as the interval on x- axis in the first equation (4), it is the time that the particle being 
in the air. The motion on the plain is equal 


t h ~ 


v 


y o 
g 


( 9 ) 


Projectile motion is generally specified the initial velocity at t = 0 isv 0 , being an angle with 
x-axis, the angle 6 0 so the components of the velocity on x-axis and y-axis is 


v *o= v o cos0 o 
v^ 0 = v o sin0 o 


( 10 ) 


Different quantities related with projectile motion which are in term of v x0 and v y() from 
equation (2) to equation (9) can be written in form of v 0 and 6 0 for instance, the range of R in 
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equation (6) in form of v 0 and 6 {) is (Bueche and Hecht, 2006; Marion and Thornton, 1995) in Fig. 
(2), Fig. (3), Fig. (4), Fig. (5). 


2 v o COS ( 6 'o) sin ( 6 l o) ' 

g 

R V o sin ( 2 ^o) 

8 


( 11 ) 


It indicates that v 0 is steady value, 6 0 angle that the highest range value is jr/4 rad. at the 
position projectile motion, with the velocity v, being angle with x-axis and y-ax is is 


v, =vcos(tf) =v^ 0 
v — vsin(#) 


The value of 0 angle is specified to be plus, being measured from +x -axis parallel to v on 
anticlockwise direction and to be minus on the opposite direction. The mentioned specification 
conforming to the direction of v y is minus. The angle value and the velocity component of projectile 

motion at different positions are indicated as the figure (2). 



Fig. 2 The velocity component v x0 and v y of projectile motion at different positions of projectile motion. 


The remark of figure (2), the velocity of x-axis is unchangeable. The size of velocity 
depends on the highest velocity of y -axis. The least of velocity size value is equal the x -axis size 
of velocity. Besides, the curve type, motion, the size of velocity and the size of angle 6 angle are 
symmetrically around the highest interval of motion. 

For example the first, numerical evaluation of the projectile motion (Bauman, 1996; 
Bauman, 2005). 

• In[l]:= eql= {m x zz [t] == 0, m y zz [t] == -mg};t == tl; 

• In[2]:= initial 1 ={x[0] == 0, x[tl] == xl, y[0] == 0, y[tl] —yl}; 

• In[3]:= dsoll = DSolve[Join[eql, initial 1 ], {x[t], y[t]}, t] // Flatten 
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—g t 2 t\ + g t t\ 2 + 2 ty\ 

2d 


_ ro , , r , tx 1 r , -gt 2 1\ +g 11\ 2 + 2ty\. 

Out[3]= {x[t] J| —, y[t] j|-2--M 

In[4]:= Map[Collect[#, t] &, dsoll, {2}] 

t x 1 st 2 As d 2 +2 yl) 

Out[4]= {x[t] J| — ,y[t] |-f- + -i-—- ’-} 

t 2 2 tl 

In[5]:= Map[Collect[#, t] &, dsoll, {2}] 

t xl 

0ut[5]= {collect[x[t], t] JJ collect[ —, t], collect[y[t], t] Jj collect[ 
t]} 

In[6]:= valuesl = {xl Jj 350, yl JJ 0, g Jj 9.8}; 

In[7]:= coordl[t_, tl] ={x[t], y[t]} /. Dsoll /. valuesl 
_ r „ , 350 t -9.8 t 2 /l + 9.8 11\ 2 

0u,[7]=< —•- 17 ^- } 

In[8]:= points[tl_]:= Table[coordl[t, tl], {t, 0, tl, 0.5}] 

In[9]:= plotl[tl_]:= ListPlot[points[tl] // Evaluate, PlotStyle JJ {Hue[0.9], PointSize[0.02]}, 
GridLines JJ Automatic, DisplayFunction JJ Identity] 

In[10]:= Show[{plotl[5], plotl[10], plotl[15], plotl[20]}, Epilog JJ {Text[“tl = 5”, 
coordl[5/2, 6.5], TextStyle JJ {FontWeight | “Bold”}], Text[“tl = 10”, coordl[10/2, 
11.5], TextStyle JJ {FontWeight JJ “Bold”}], Text[“tl = 15”, coordl[15/2, 16.0], TextStyle JJ 
{FontWeight JJ “Bold”}], Text[“tl = 20”, coordl[20/2, 19.5], TextStyle ]) {FontWeight JJ 
“Bold”} ] }, DisplayFunction J| $ DisplayFunction, PlotRange | All] 

Out[10]= graph see Fig. (3), Fig. (4) and Fig. (5) 
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Fig. 3 A projectile(a ball) moving under the force of gravity is at its maximum height when t\ = 20 and 
V = 0 and minimum height when t\ = 5 . At that moment, the ball is traveling horizontally. 
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Fig. 4 A projectile(a ball) moving under the force of gravity is at its maximum height when tl = 20 s. and 
vy — 0 m/s and minimum height when t\ = 5 s. . By Setting xl =350 m., y 1=100 m., g=9.8m/s 2 . 



Fig. 5. A projectile(a ball) moving under the force of gravity is at its maximum height when tl = 20s. and 
Vy = 0 m/s and minimum height when tl = 5 s.. By Setting xl =350m, y1=200m, g=9.8m/s 2 . 

For example the second, numerical evaluation of the projectile motion (Wolfram, S. 
1994; Maeder, R. 1991). 

• In[ll]:= initial = {x[0] == 0, x z [0] == vO Cos[0], y[0] == 0, y z [0] == vO Sin[tf ]}; 

• In[12]:= dsol2 = DSolve[Join[eql, initial2], {x,y} , t] // Flatten 

1 9 

• Out[12]= {x JJ Function[{t} , t vO Cos[#]] , y Jj Function[{t} , — (-g t 2 + 2 t vO Sin[<9])]} 

• In [13]:= eq2 = y z [t] == 0/. Dsol2 
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Out[13]=| (- 2 g t + 2 vO Sin[0]) == 0 

In[14]:= y z [t] == 0/. Dsoll 

Out[14]= y z [t] == 0 

In[15]:= tsol = Solve[eq2, t] // Flatten 

Ou,[15]= ( t^^«} 

g 


In[16]:= {x[t] , y[t]} /. Dsol2 /. tsol // Simplify 


Out[16]= { 


v0 2 Cos[<9]Sin[6>] v0 2 (Sin[6>]) 2 


g 


} 


In[17]:= values2 = {vO J| 100, g JJ 9.8}; 

In[18]:= coord2[t_, 6 ={x[t], y[t]} /. Dsol2 /. Values2 


Out[18]= {100 t Cos[0],- (-9.8 t 2 + 200 t Sin[0])} 

In[19]:= plot2[#_]:= ListPlot[Table[coord2[t, d ~\, {t, 0,20 , /4}] // Evaluate, PlotStyle JJ 
{Hue[0.9], PointSize[0.02]}, GridLines JJ Automatic, DisplayFunction JJ Identity] 

In[20]:= Show[{plot2[^], plot2[^J, plot2[y ], plot2[^]}, Epilog J| {Text[“0 = 

3 6 4 12 3 


coord2[4.5, ^-], TextStyle JJ {FontWeight JJ “Bold”}], 


Text[“# = — ”, coord2[5, — ], 
6 6 


TextStyle JJ {FontWeight JJ “Bold”}], Text[“# = coord2[5, 7 E\, TextStyle JJ 

{FontWeight JJ “Bold”}], Text[“# = coord2[5, TextStyle JJ {FontWeight JJ 

“Bold”} ] }, DisplayFunction JJ $ DisplayFunction, PlotRange JJ {{0 , 1100} , {0,500}} , 
FrameLabel J| {“x[t]” , “y[t]”}] 

• Out[20]= graph see figure (4) (Ginsberg, E.S. 1995) 


y 



Fig. 6. For a fixed initial velocity (ball) and if air resistance is ignored, a projectile's trajectory will have a 
maximum range for an elevation angle of jtj 4 rad . (Ginsberg, E.S. 1995). The range is the horizontal 
distance the projectile travels to reach the same height from which it started. 
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Fig. 7. For a fixed initial velocity (ball) and if air resistance is ignored , a projectile’s trajectory will have a 
maximum range for an elevation angle of jt/ 4 rad. (Ginsberg , E.S. 1995). From this is Fig. 7, if the value 
initial velocity (VO) has increase, the range is the horizontal distance the projectile travels to has lessen. By 

setting V0=90 m/s. 
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Fig. 8. For a fixed initial velocity (ball) and if air resistance is ignored , a projectile’s trajectory will have a 
maximum range for an elevation angle of jt/ 4 rad. (Ginsberg, E.S. 1995). From this is Fig. 8, if the value 
initial velocity (VO) has increase, the range is the horizontal distance the projectile travels to has lessen. By 

setting V0=80 m/s 


For example the third, numerical evaluation of the projectile motion (Wolfram, 1994; 
Maeder, 1991). 
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ln[l]:= Manipulate [ParametricPlot [ { {v Cos [0] t, VSin[0] t-9.8t 2 /2}}, 

{t, 0, tf}, PlotStyle {{Thick, Black}}, AxesOrigin {0, 0}, 

FrameLabel^ {Style[x. Bold, Black, Large], Style[y. Bold, Black, Large]}, 
LabelStyle Directive[Black, Bold], LabelStyleDirective[Black, Bold], 
Frame -► True, Frame Style -¥ Directive [Thick, Black, 16] , Axes True, 
PlotRange -* {{0, 400}, {0, 300}}, ImageSize 500, Epilog -* 

{Black, PointSi ze [.02] , Point [ { {v Cos [0] tf, V Sin [0] tf-9.8tf 2 /2jj]J], 
{{V, 48.8, "initial velocity"}, 1, 100, Appearance 4 "Labeled"}, 

{{0, 1.09, "angle"}, .1, tt/ 2. Appearance "Labeled"}, 

{{vt, 100, "terminal velocity"}, .01, 300, Appearance # "Labeled"}, 

{{tf, 8.75, "time"}, .1, 17, Appearance"Labeled"}] 
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Fig. 9: A projectile fire from the origin with an initial speed of 60 m/s at various time of 9 second. 
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Method 

The researcher divided the samples into 2 groups, experiment group and controlled group by 
random. Fifty students of experiment group study projectile motion by using Mathematica program 
and fifty students of controlled group study projectile motion by using projectile motion 
experiment. 

Get both experiment group and controlled group to take pre-test by using 30 items of 
achievement test. 

Let experiment group study projectile motion, using Mathematica program and controlled 
group study projectile motion, using projectile motion experiment set. 

Get both experiment group and controlled group to take 60 minutes post-test by using 30 
items of achievement test. 

Instruments 

The instruments used for projectile motion experiment are Mathematica program and projectile 
motion experiment set. 

The instrument used for collecting data is the 30 items of achievement test. 

Data analysis 

Analyze the physics achievement of Faculty of Science and Technology Phetchabun Rajabhat 
University students studying projectile motion by using Mathematica program and Faculty of 
Science and Technology Phetchabun Rajabhat University students, studying projectile motion by 

using experiment set of projectile motion, finding the average (x) (The simplest number used to 
characterize a sample is the mean, which for N values x i , i = 1,2,K , N.) is defined by (Riley, K.F. 
et al., 2006) 



(13) 


and the sample standard deviation (SD) is the positive square root of the sample variance, i.e. 



We may therefore write the sample variance SD 2 as 



from which the sample standard deviation is found by taking the positive square root. Thur, by 
evaluating the quantities and ^ ^ xf for our sample, we can calculate the sample mean 

and sample standard deviation at the same time. 
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Compare pre-test and post-test (Lord, F.M. 1967; Airasian, P.W. and Madaus, G.F., 1972) of 
physics achievement of Faculty of Science and Technology Phetchabun Rajabhat University 
students, studying projectile motion experiment by using Mathematica program and Faculty of 
Science and Technology Phetchabun Rajabhat University students studying projectile motion by 
using and experiment set of projectile by t-test independent. 

The Results 

The results of analyzing Pre-test and post-test physics achievement; projectile of Faculty of Science 
and Technology Phetchabun Rajabhat University students are as the following. From table (1) 
finding an average of the pre-test of physics achievement, projectile motion of experiment group, 
studying projectile motion experiment by using Mathematica program is 12.04 and an average of 
the pre-test of physics of controlled group studying projectile motion by using projectile motion 
experiment set is 13.30. Using t-test independent Sample to compare achievement finds that the 
experiment group and controlled group’s achievement is different insignificantly. From table (2), 
Finding the average of physics achievement, projectile motion by using Mathematica program is 
24.38 and the average of controlled group studying projectile motion by using projectile motion 
experiment set is 20.82 Using t-test independent sample to compare achievement finds that 
controlled group studying projectile motion by using Mathematica program is significantly higher at 
level of 0.01. It means the achievement of Faculty of Science and Technology Phetchabun Rajabhat 
University students studying projectile motion by using Mathematica program is higher than the 
achievement of Faculty of Science and Technology Phetchabun Rajabhat University students 
studying projectile motion by using projectile motion. 


Table 1. The result of physics achievement study, projectile pre-test 


Statistics 

Sample 

N 

X 

SD 

t-test 

Experiment group 

50 

12.04 

2.821 

3.743 

Controlled group 

50 

13.30 

1.982 



Table 2. The result of physics achievement study, projectile post-test 


Statistics 

Sample 

N 

X 

SD 

t-test 

Experiment group 

50 

24.38 

1.772 

12.04 

Controlled group 

50 

20.82 

1.289 



Conclusion and Discussion 

The achievement of Faculty of Science and Technology Phetchabun Rajabhat University students 
studying projectile motion by using Mathematica is higher than the achievement of Faculty of 
Science and Technology Phetchabun Rajabhat University students studying projectile motion by 
using projectile motion experiment set at level of 0.01. 

The achievement of Faculty of Science and Technology Phetchabun Rajabhat University 
Technology students who use Mathematica program is higher than the achievement of Faculty of 
Science and Technology Phetchabun Rajabhat University students studying projectile motion by 
using projectile motion experiment set. Examining hypothesis finds that the achievement of Faculty 
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of Science and Technology Phetchabun Rajabhat University Technology students using 
Mathematica program is higher than the achievement of Faculty of Science and Technology 
Phetchabun Rajabhat University students studying projectile motion by using projectile motion 
experiment set at level of 0.01 which relates to the specified hypothesis of study since Mathematica 
program is an easy program. It can be experimented at home or at university. It inspires students to 
study because it’s strange and new for studying which related to the act of education 1999 B.E. part 
4 section 22-23 to encourage students to study naturally and efficiently, let students study 
autonomously, encourage students to enjoy learning, to understand lessons well and make students 
creative as the purposes of study. 
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